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ABSTRACT 

We present the results of the numerical calculation of the probability distribution of 
the value of the monopole creation operator in SU(2) lattice gluodynamics. We work 
in the maximal abelian projection. It occurs that at the low temperature, below 
the deconfinement phase transition the maximum of the distribution is shifted from 
zero, which means that the effective constraint potential is of the Higgs type. Above 
the phase transition the minimum of the potential (the maximum of the monopole 
field distribution) is at the zero value of the monopole field. This is the direct proof 
of the existence of the abelian monopole condensate in the confinement phase of the 
gluodynamics, which confirms the dual superconductor model of the confining vacuum. 



1. Introduction 

The monopole mechanism of the colour confinement is generally accepted by the commu- 
nity. Still there are many open questions. In the lattice gluodynamics it is very important 
to find the order parameter, constructed from the monopole field, for the deconfinement 
phase transition. The first candidate is the monopole condensate, which should be nonzero 
in the confinement phase and vanish at the phase transition. To study the monopole con- 
densate we need the explicit expression for the operator & m0 n(%), which creates the abelian 
monopole at the point x. The operator &mon{x) was found for the compact electrodynamics 
with the Villain form of the action by Frohlich and Marchetti Q , and it was studied numer- 
ically in refs.0. In Section 2 we construct the monopole creation operator for an arbitrary 
abelian projection of lattice SU(2) gluodynamics. The numerical results presented in Sec- 
tion 3 are obtained for the maximal abelian projection, for this projection many numerical 
simulations show that the gluodynamic vacuum behaves as the dual superconductor (see 
H and references therein). In refs.|Q] the another form of the monopole creation operator 
was studied, and it was found that its expectation value vanishes in the deconfinement 
phase; as we discuss at the end of Section 2, our operator is positively defined, therefore 
our definition of ^mon(x) differs from that of ref.Q, still our conclusions and that of ref.|Q] 

* Talk given at 29-th International Symposium on the Theory of Elementary Particles, Buckov, 29 August 
- 2 September, 1995. 



coincide: the monopole condensate exists in the confinement phase of lattice gluodynamics. 
The analogous claim is done in ref.||], where the monopole condensate is calculated on the 
basis of the percolation properties of the monopole currents. 



2. Monopole Creation Operator 

First we give the formal construction of the monopole creation operator for the abelian 
projection of SU{2) gluodynamics. Let us parametrize SU(2) link matrix in the standard 
way: = cos^e^M; U™ = sin0^e^; u™ = U% = -U%*; < < 

ir/2, -tt < 9,x < 7T. 

The plaquette action in terms of the angles 0, and \ can De written as follows: 

S P = ^Tr U 1 U 2 U+Ut = S a + S n + S\ (1) 

where 

S a = COsOp COS(/)l COS 02 COS 03 COS 04, (2) 

S n and S l describe the interaction of the fields 6 and x an d selfinteraction of the field x 0> 

dp = 61 + 9 2 — 03 — 6*4 , (3) 

here the subscripts 1, 4 correspond to the links of the plaquette: 1 — > {x, x + /}}, 4 — > 
{x, x + z>}. 

If we fix the abelian projection, each term S a , S n and S l is invariant under the residual 
U(l) gauge transformations: 



+ a x - a x+ fx , (4) 
Xxfi -> Xxfi + a x + a x+ f t . (5) 

The operator which creates the monopoles at the point x of the dual lattice is defined 
as follows: 



U{x) = exp {P[-S(9 P , 0) + S(0 P + W P (x), 0)]} , (6) 
we define the function Wp(x) below. Substituting eq.(|l|-[|) in eq.(||) we get 



U(x) = exp|^/3[-cos(0p) + cos(0p + Wp(x))]| , 



(7) 



where (3 = cos 0i cos 02 cos 03 cos 04 (3. Effectively the monopole creation operator shifts all 
abelian plaquette angles 0p. 



2 



For the compact electrodynamics with the Villain type of the action the above definition 
coincides with the definition of Frdhlich and Marchetti ||]. For the general type of the action 
in compact electrodynamics we can use the described above construction. The proof is given 
in Appendix A. The gluodynamics in the abelian projection contains the compact gauge 
field 6 and the charged vector field x- The action (|l]) in terms of the fields 9 and \ is 
rather nontrivial, and at the moment we have no proof that the above construction of the 
monopole creation operator is valid in this case. Still for the rather similar Abelian - Higgs 
model, with the general type of the action, the proof exists, and it is analogous to one given 
in Appendix A. Moreover the numerical results, presented in the next section, clearly show 
that the suggested operator is the order parameter for the deconfinement phase transition. 



3. Numerical Results 

We present the results of the numerical calculations on the lattice 10 3 • 4, we impose 
the anti-periodic boundary conditions in space directions, since the construction of the 
operator IA can be done only in the time slice with the anti-periodic boundary conditions. 
The periodic boundary conditions are forbidden due to the Gauss law; formally there is no 



solution of equation (15) in the finite box with the periodic boundary conditions. To see that 
we have the order parameter for the deconfinement phase transition it is very convenient 
to study the probability distribution of the operator U. It means that we calculate the 
expectation value < 5(ip — U{x)) >. The quantity like the effective constraint potential, 



-^)=<%-l£W(z))> (8) 



has more physical meaning than the probability distribution. At the moment we have no 
enough statistics to calculate V e ff(ip), and we present our results for the quantity V(ip), 
defined as: 



e -v( v ) =< s{ip-U{x)) >. (9) 

In Figs. 1(a) and 1(b) we show V{ip) for the confinement and the deconfinement phases. 
It is clearly seen that in the confinement phase the minimum of V(<p) is shifted from zero, 
while in the deconfinement phase the minimum is at zero value of the monopole field (p. We 
used the positively definite operator U(x) @, but in the dual representation the creation 
operator of the monopole ( |l~4|) is not positively defined, the sign is loosed since we perform 
the inverse duality transformation on the infinite lattice. On the finite lattice it is possible 
to get the non-positive defined operator U(x), in that case instead of Fig. 1(a) we get the 
Higgs - type potential. This little bit more complicated calculations are now in progress. 
Still Figs. l(a),(b) clearly show that the position of the minimum of V(tp) plays the role of 
the order parameter. On Fig. 2 we show the dependence of the position of the minimum, 
<p c , on the temperature, it is seen that ip c vanishes at the point of the phase transition. 
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Appendix A 

Below we construct the monopole creation operator for the compact electrodynamics 
with the general type of the action; the similar construction exists for the compact Abelian- 
Higgs model with the general type of the action. First we perform the duality transformation 
of the partition function for the AD lattice compact electrodynamics: 

Z = J V9exp{-S(d9)}, (10) 

— 7T 

We use the notations of the calculus of differential forms on the lattice J/]] (see also 
Appendix B). The symbol J T>9 denotes the integral over all link variables 9. The partition 
function of the dual theory, 



Z d = £ exp{-*S(d*n)}, (11) 

*n(c 1 )e^ 

+7T 

*S(p) = -In J VFexp{-S(F)+i(F,p)}, (12) 

— 7T 

can be represented as the following limit of the partition function for the Abelian-Higgs 
theory: 



+TT +00 

Z d = Jim^J V*(p J V*B J2 ex.p{-*S(d*B/2ir) - n\\*B -d*ip + 2ir*n\\ 2 }, (13) 

— 7v -oo *n(ci)e% 

here *S(d*-B/27r) is the kinetic energy of the dual gauge field *B (the analogue of F^ u ) 
and the Higgs field exp{i*ip} carry magnetic charge, since it interacts via the covariant 
derivative with the dual gauge field *B. The Dirac operator H, 



U d (x) = e i ^-exp{-i(*D x ,*B)}, (14) 

5*D X = *5 X (15) 
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is the gauge invariant monopole creation operator. It creates the cloud of photons and the 



monopole at the point x. In (15) *5 X is the lattice 5-function, it equals to unity at the cite 
x of the dual lattice and is zero at other cites. Note that in the above formulas the radial 
part of the Higgs field which carry the magnetic charge is fixed to unity. 

Coming back to the original partition function ( JTo| ) we get the expectation value of the 
monopole creation operator in terms of the fields 9: 

<U{x)>= — J Vdexp{-S(de + W P (x))}, (16) 

— 7T 

W P (x) = 2tt5A-\D x - w x )), (17) 

where the Dirac string attached to the monopole is represented by the integer valued 
1-form *u> x , which satisfies the equation: 5*uj x = *5 X . 

The partition function ( [Tot) for the electrodynamics with the general type of the action 
can be represented as the sum over the monopole closed currents Q. It is straightforward to 
show that the monopole creation operator U{x) creates the not-closed monopole trajectory 
which starts at the point x. This fact shows that lA{x) is the monopole creation operator. 



Appendix B 

Below we briefly summarize the main notions from the theory of differential forms on 
the lattice [fjj. The advantages of the calculus of differential forms consists in the general 
character of the expressions obtained. Most of the transformations depend neither on 
the space-time dimension, nor on the rank of the fields. With minor modifications, the 
transformations are valid for lattices of any form (triangular, hypercubic, random, etc). A 
differential form of rank k on the lattice is a function defined on /c-dimensional cells 
Cfc of the lattice, eg the scalar (gauge) field is a 0-form (1-form). The exterior differential 
operator d is defined as follows: 



(d$( Cfc+1 )= Yl ( 18 ) 

Here dck is the oriented boundary of the £>cell c&. Thus the operator d increases the rank of 
the form by unity; dip is the link variable constructed, as usual, in terms of the site angles p, 
and dA is the plaquette variable constructed from the link variables A. The scalar product 
is defined in the standard way: if ip and i\) are fc-forms, then (ip,ip) = J2c k L P{ c k)'4 ) {.Ck)^ 
where J2c k ^ s the sum over all cells c^. To any /c-form on the D-dimensional lattice there 
corresponds a [D — /c)-form *&(*Ck) on the dual lattice, *Ck being the (D — k) -dimensional 
cell on the dual lattice. The codifferential 5 = *d* satisfies the partial integration rule: 
(p, 5ip) = (dip, tp). Note that 5&(ck) is a (k — l)-form and 5<&(cq) = 0. The norm is defined 
by: ||a|| 2 = (a, a); therefore, \\B — dp + 27m|| 2 in (13) implies summation over all links. 
J2i(ci)e& denotes the sum over all configurations of the integers / attached to the links c\. 
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The action (|T^) is invariant under the gauge transformations B' = B + da, ip' = 99 + a due 
to the well known property d 2 = 5 2 = 0. The lattice Laplacian is defined by: A = d5 + 5d. 



Figure Captions 

Figl. V{ip) for the confinement (a) and the deconfinement (b) phases. 
Fig2. Position of the minimum ip c of V(ip) vs. temperature T. 
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